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Excitations of three-dimensional spin glasses are computed 
numerically. We find that one can flip a finite fraction of an 
L x L x L lattice with an O(f) energy cost, confirming the 
mean field picture of a non-trivial spin overlap distribution 
P(q). These low energy excitations are not domain- wall-like, 
rather they are topologically non-trivial and they reach out to 
the boundaries of the lattice. Their surface to volume ratios 
decrease as L increases and may asymptotically go to zero. If 
so, link and window overlaps between the ground state and 
these excited states become "trivial" . 
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Spin glasses jlj] are currently at the center of a hot de- 
bate. One outstanding question is whether there exists 
macroscopically different valleys whose contributions si- 
multaneously dominate the partition function. At zero 
temperature, given the ground state configuration, this 
leads one to ask whether it is possible to flip a finite frac- 
tion of the spins and reach a state with excess energy 
0(1). From a mean- field perspective one expects this 
to be true since it happens in the Sherrington-Kirkpatrick 
(SK) model. However it is very unnatural in the con- 
text of the droplet || or scaling [Q approaches where 
the characteristic energy of an excitation grows with its 
size. Recently it has been suggested that the energy of 
an excitation may grow with its size £ as in the droplet 
scaling law, E{£) w £ 01 , but only for f <C L, and that 
for £ = 0(L) the energies cross over to a different law, 
E ps L 8g , where L is the size of the system j5|. The 
first exponent, 9i (I for local), may be given by domain 
wall estimates, Oi 0.2, while the second exponent, 6 g (g 
for global), could be given by the mean field prediction, 
g = 0. In this "mixed" scenario, one has coexistence of 
the droplet model at finite length scales and a mean-field 
behavior (if 9 g = 0) for system-size excitations (£ w L for 
which a finite fraction of all the spins are flipped) . 

The purpose of this article is to provide numerical evi- 
dence that such a mixed scenario is at work in the three- 
dimensional Edwards- Anderson spin glass. We have de- 
termined ground states and excited states for different 
lattice sizes and have analyzed their geometrical prop- 
erties. The qualitative picture we reach is that indeed 
g w 0. System-size constant energy excitations are 
not artefacts of trapped domain walls caused by periodic 
boundary conditions, they are intrinsic to this kind of 
frustrated system. The energy landscape of the Edwards- 
Anderson model then consists of many valleys, probably 



separated by large energy barriers. Extrapolating to fi- 
nite temperature, this picture leads to a non-trivial equi- 
librium spin overlap distribution function P(q). 

Given the geometric properties of our excitations, we 
suggest a new scenario for finite dimensional spin glasses: 
if the surface to volume ratios of these large scale exci- 
tations go to zero in the large L limit, then the replica 
symmetry breaking will be associated with a trivial link 
overlap distribution function P(qi). We have coined this 
scenario TNT for trivial link overlaps yet non-trivial spin 
overlaps. Such a departure from the standard mean field 
picture might hold in any dimension d > 3. 

The spin glass model - We consider an Edwards- 
Anderson Hamiltonian on a three-dimensional L x L x L 
cubic lattice: 
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The sum is over all nearest neighbor spins of the lat- 
tice. The quenched couplings Jy are independent ran- 
dom variables, taken from a Gaussian distribution of zero 
mean and unit variance. For the boundaries, we have 
imposed either periodic or free boundary conditions. Al- 
though in simulations of most systems it is best to use 
periodic boundary conditions so as to minimize finite size 
corrections, the interpretation of our data is simpler for 
free boundary conditions. It may also be useful to note 
that if boundary conditions matter in the infinite volume 
limit, free boundary conditions are the experimentally 
appropriate ones to use. 

Extracting excited states - The problem of finding 
the ground state of a spin glass is a difficult one. In this 
study we use a previously tested Q algorithmic proce- 
dure which, given enough computational ressources, gives 
the ground state with a very high probability for lattice 
sizes up to 12 x 12 x 12. (Since our J,jS are continuous, 
the ground state is unique up to a global spin flip.) Our 
study here is limited to sizes L < 11; then the rare er- 
rors in obtaining the ground states are far less important 
than our statistical errors or than the uncertainties in 
extrapolating our results to the L — -> oo limit. 

Our purpose is to extract low-lying excited states to 
see whether there are valleys as in the mean field picture 
or whether the characteristic energies of the lowest-lying 
large scale excitations grow with L as expected in the 
droplet/scaling picture. Ideally, one would like to have a 
list of all the states whose excess energy is below a given 
cut-off. However, because there is a non-zero density 
of states associated with droplets (localized excitations) , 
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this is an impossible task for the sizes of interest to us. 
Thus instead we extract our excitations as follows. Given 
the ground state (hereafter called Co), we choose two 
spins Si and Sj at random and force their relative orien- 
tation to be opposite from what it is in the ground state. 
This constraint can be implemented by replacing the two 
spins by one new spin giving the orientation of the first 
spin, the other one being its "slave". We then solve for 
the ground state C of this modified spin glass. The new 
state C is necessarily distinct from Co as at least one spin 
(Si or Sj) is flipped. That flipped spin may drag along 
with it some of its surrounding spins, forming a droplet 
of characteristic energy 0(1). In the droplet picture, this 
is all that happens in the infinite volume limit. However 
if there exist large scale excitations with 0(1) energies, 
then C may be such an excitation if its energy is below 
that of all the droplets containing cither Si or Sj . 

Statistics of cluster sizes - Let V be the number 
of sites of the cluster defining the spins that are flipped 
when going from Cq to C (by symmetry, V is taken in 
[1, L 3 /2)). If P(V) is the probability to have an event of 
size V, the droplet and mean field pictures lead us to the 
following parametrization: 

P(V) = (l-a)P l (V)+aP g (V/L 3 ). (2) 

Here, Pi and P g are normalized probability distributions 
associated with the droplet events (V fixed, L — > oo) and 
the global events (V = 0(L 3 )). If large scale excitations 
have energies 0(L° a ), the ratio a /(I — a) of the two con- 
tributions should go as L~ e «. In the droplet/scaling pic- 
ture, the global part decreases as L~ Sl ; that is slow since 
6i w 0.2. In contrast, in the mean field scenario, both the 
V finite and the V growing linearly with L 3 contributions 
converge with non-zero weights, < a < 1, albeit with 
0(L~ 6 ') finite size corrections. 

Given that the usable range in L is no more than a 
factor of two so that L~ Bl does not vary much, measure- 
ments of P(V) on their own are unlikely to provide strin- 
gent tests. Nevertheless, consider the probability Q(v, v') 
that V/L 3 is in the interval [«,«']. Up to finite size cor- 
rections, Q(v,v') = a J P g (x)dx. In our computations, 
we have used 5 < L < 11, averaging for each L over 
2000 to 10000 randomly generated samples of the Jij. 
For each sample, we determined the ground state, and 
then obtained 3 excitations by choosing successively at 
random 3 pairs of spins (Si,Sj). We find that Q(v, v') de- 
creases slowly with L for both periodic and free boundary 
conditions, as expected in the droplet and mean field pic- 
tures. Because 9i is small, when we perform fits of the 
form Q(v,v') = A + BL~^ : we are not able to exclude 
A = nor A ^ with any significant confidence, so a 
more refined method of analysis is necessary: we will thus 
consider the geometrical properties of the events. 

Before doing so, note that the statistical error on 
Q(v,v') depends on the number of large scale events 




FIG. 1. Example of excitation found for a 12 lattice with 
free boundary conditions. 

found in the [v,v'] interval. If the spin Si or Sj has a 
small local field, there is a good chance that the corre- 
sponding event will have V = 1, thereby reducing the 
statistics of the "interesting" events. To amplify our sig- 
nal of large V events, we did not consider such spins and 
focused our attention on spins in the top 25 percentile 
when ranked according to their local field. All of our 
data was obtained with that way of selecting Si and Sj . 
(Naturally, P(V) and Q(v, v') depend on this choice, but 
the general behavior should be the same for any choice.) 

Topological features of the clusters — Our claim that 
9 g < 8i can be credible only if our large scale excita- 
tions arc different from domain-walls (whose energies are 
believed to grow as L - 2 ). It is thus useful to consider ge- 
ometrical characterizations of the excitations generated 
by our procedure. Figure |l| shows a typical cluster found 
for a 12 3 lattice. It contains 622 spins and its (excitation) 
energy is 0.98 which is O(l). The example displayed is 
for free boundary conditions which permits a better vi- 
sualization than periodic boundary conditions. 

The cluster shown touches many of the 6 faces of the 
cube, and the same is true for the complement of that 
cluster. Such a cluster has a very non-trivial topology 
and is thus very far from being domain- wall like. This 
motivates the following three-fold classification of the 
events we obtain when considering free boundary con- 
ditions. In the first class, a cluster and its complement 
touch all 6 faces of the cube. In the second class, a clus- 
ter touches at most 3 faces of the cube. The third class 
consists of all other events. Finite size droplets should 
asymptotically always fall into the second class, albeit 
with finite size corrections of order L~° l . 

Does the first class constitute a non-zero fraction of 
all events? At finite L, we find the following fractions: 
23.3% (L = 5), 23.9% (L — 6), 25.1% (L — 7), 24.4% 
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FIG. 2. Integrated probability Q(v, 1/2) of events in the 
first class. (From bottom to top, L = 5, 7, 9, 11. 

(L = 8), 25.0% (L — 9), 25.7% (L = 10), and 26.0% 
(L = 11). The trend of these numbers suggests that 
the first class does indeed encompass a finite fraction of 
all the events when L — > oo. We also considered the 
scaling of cluster sizes with L. Fig. || shows Q(v, 1/2) as 
a function of « = ^/i 3 , restricted to events belonging 
to the first class. (The v — values are the fractions 
we just gave above.) The curves for different L show 
a small drift, Q(v, 1/2) growing with L. We consider 
this drift to be a finite size effect and that the correct 
interpretation of our data is 9 g ps 0, in agreement with 
the mean field picture. Our conclusion is then that as 
L —> oo, there is a finite probability of having an 0(1) 
energy excitation that is non-domain- wall like, the cluster 
and its complement touching all faces of the cube. 

Surface to volume ratios - Obviously the mean field 
picture obtained by extrapolating results from the SK 
or Viana-Bray spin glasses cannot teach us anything 
about the topology of excitations for three-dimensional 
lattices. But mean field may serve as a guide for 
other properties such as the link overlap qi between 
ground states and excited states. In the SK model, 
the spin overlap q = ^2 SiS'jN and the link overlap 
91 = E(S l S J )(S' l S' :j )/(N(N - l)/2) satisfy Ql = q 2 , and 
both q and qi have non-trivial distributions. Extrapo- 
lating this to our three dimensional system, the mean 
field picture predicts that the clusters associated with 
large scale excitations both span the whole system (as 
we saw with free boundary conditions) and are space fill- 
ing. Quantitatively, this implies that their surface grows 
as the total volume of the system, i.e., as L 3 . 

To investigate this question, we have measured the 
surface of our excitations, defined as the number S of 
links connecting the corresponding cluster to its complc- 
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TABLE I. Chi squared per degree of freedom for the fits to 
the data of Fig. |[ 



ment. (Then qi = 1 — 2S/3L 3 .) In figure || we show 
the mean value of S/L 3 as a function of L for v — V/L 3 
belonging to the three intervals ]0.20,0.25], ]0.30,0.35], 
and ]0. 40, 0.45]. The data shown are for free boundary 
conditions, but the results are very similar for periodic 
boundary conditions. The most striking feature is that 
the curves decrease very clearly with L. For each inter- 
val, we have fitted the data to (S)/L 3 = A + B/L^ and 
to a polynomial in 1/L. Of major interest is the value 
of the constant because it gives the large L limit of the 
curves. 

Table | summarizes the quality of the fits as given by 
their Xr ( cn i squared per degree of freedom) . In all cases 
the fits are reasonably good; this is not so surprizing 
because our range of L values is small. The most reliable 
fits are obtained using a quadratic polynomial in 1/L, 
this functional form leading to a smooth and monotone 
behavior of the parameters and to small uncertainties in 
the parameters. For the large L limits, these fits give 
A = 0.22, A = 0.27 and A = 0.30 for the three intervals. 
(We do not give results for the linear fits which on the 
contrary are very poor.) The constant plus power fits 
also have good Xr but the As obtained were small and 
decreased with v; also they had large uncertainties and 
seemed to be compatible with A — 0. Because of this, we 
also performed fits of the form (S) / L 3 = B /L^ . These 
are displayed in Fig. ||and lead to \i « 0.30 (the exponent 
varies little from curve to curve), again with reasonable 
X^s. Because of this, we feel we cannot conclude from 
the data that the surface to volume ratios tend towards 
a non-zero asymptote. What can be said is that this 
asymptote seems to be small, and that it will be difficult 
to be sure that it is non-zero without going to larger 
values of L. 

Discussion - For the three dimensional Edwards- 
Anderson spin glass model, we have presented numeri- 
cal evidence that it is possible to flip a finite fraction of 
the whole lattice at an energy cost of O(l), correspond- 
ing to 6 g rs as predicted by mean field. This property 
transpired most clearly through the use of free boundary 
conditions, allowing one to conclude that 6 g « is not 
an artefact of trapped domain walls caused by periodic 
boundary conditions. Extrapolating to finite tempera- 
ture, we expect the equilibrium P(q) to be non trivial as 
in the mean field picture. 

The other messages of our work concern the nature 
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FIG. 3. Mean value of surface to L 3 ratios for v = V/L 3 
in intervals around 0.225, 0.325, 0.425 (bottom to top) using 
free boundary conditions. Curves are pure power fits. 



of these large scale excitations whose energies are 0(1). 
First, using free boundary conditions, we found them to 
be topologically highly non-trivial: with a finite prob- 
ability they reach the boundaries on all 6 faces of the 
cube. Thus they are not domain-wall-like, rather they 
are sponge-like. Second, our data (both for periodic and 
free boundary conditions) indicate very clearly that their 
surface to volume ratios decrease as L increases. The 
most important issue here is whether or not these ra- 
tios decrease to zero in the large L limit. Although our 
data are compatible with a non-zero limiting value as 
predicted by mean field, the fits were not conclusive so 
further work is necessary. 

If the surface to volume ratios turned out to go to zero, 
we would be lead to a new scenario that we have coined 
"TNT". In the standard mean field picture, the surface 
to volume ratios cannot go to zero; indeed in the SK and 
Viana-Bray spin glass models there are no spin clusters 
with surface to volume ratios going to zero. However, in 
finite dimensions, one can have surface to volume ratios 
going to zero, in which case qi — ► 1. This property would 
then lead to a non-trivial P(q) but to a trivial P{qi)- 
This trivial-non-trivial (TNT) scenario does not seem to 
have been proposed previously. 

Perhaps the most dramatic consequence of this new 
scenario is for window overlaps in spin glasses: because 
in TNT one is asymptotically always in the bulk of an 
excitation, correlation functions at any finite distance 
will show no effects of replica symmetry breaking. That 
this may arise in fact is supported by work by Palassini 
and Young Q who showed that certain window overlaps 
seemed to become trivial as L — > oo. (See also § for 
a similar discussion in two-dimensions.) These authors 



referred to this property as evidence for a "trivial ground 
state structure". But in our picture the global (infinite 
distance) structure is non- trivial, as indicated by 9 g = 0, 
in sharp contrast to the droplet/scaling picture. Also, 
in very recent work Q, Palassini and Young have ex- 
tended their previous investigations and have extracted 
excited states by a quite different method from ours, and 
they find that their data is compatible with the TNT 
scenario. Naturally, there is also evidence in favor of 
the non-triviality of window overlaps |u|. Nevertheless, 
we believe that our mixed scenario is a worthy candi- 
date to describe the physics of short range spin glasses. 
Furthermore, its plausibility should not restricted to 3 
dimensions, it could hold in all dimensions greater than 
2. (Note that in d = 2, excitations are necessarily topo- 
logically trivial.) An important indication of this was 
obtained by Palassini and Young whose computations Q 
favor the TNT scenario over the droplet picture in the 
4-dimensional Edwards- Anderson model. 
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